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ABSTRACT 

Finite-amplitude  standing  waves  in  air  at  ambient  tem- 
peratures contained  within  a  rigid-walled  rectangular  cavity 
having  one  variable  interior  dimension  were  experimentally 
investigated o   The  effect  of  degenerate  or  nearly  degenerate 
configurations  on  the  harmonic  content  of  the  observed  pres- 
sure waveform  was  compared  to  the  theory  of  Coppens  and 
Sanders  as  modified  to  include  empirically  determined  losses 
and  resonances.   Second  and  third  harmonic  distortion  for 
both  degenerate  and  nondegener a te  cavity  configurations  were 
recorded  as  continuous  functions  of  frequency  in  the  vicinity 
of  the  resonances  for  the  100  and  110  modes0   It  was  found 
that  in  rectangular  cavities  having  walls  free  of  any  per- 
turbation, the  theoretical  model  can  accurately  predict  both 
shapes  and  amplitudes  of  the  harmonic  content  of  the  finite- 
amplitude  standing  wave  whether  or  not  degeneracies  are 
present. 
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I .    INTRODUCTION 

In  1968,  Coppens  and  Sanders  [1]  applied  a  one-dimen- 
sional, nonlinear,  acoustic  wave  equation  with  a  dissipa- 
tive  terra  describing  absorptive  losses  to  r igid_walled , 
closed  tubes  with  large  length-to-diameter  ratios.   They 
later  expanded  this  theory  to  incorporate  empirically 
determined  losses  and  resonance  frequencies  [2] .   The  model 
was  further  extended  to  include  two-  and  three-dimensional 
cases,  which  was   experimentally  investigated  by  Lane  [3] 
and  Devall  [4 ] . 

One  of  the  significant  results  of  these  latter  investi- 
gations was  that  the  agreement  between  theory  and  experi- 
ment deteriorated  when  degenerate  modes  were  present.   The 
theory  was  not  able  to  account  for  the  effect  of  modes 
degenerate  to  members  of  the  family  of  the  driven  mode. 

To  study  the  effects  of  degeneracies,  a  rectangular 
cavity  was  designed  and  constructed  so  that  one  wall  could 
be  moved  and  secured  at  various  positions  to  introduce  or 
remove  degenerate  modes  as  desired.   The  primary  interest 
in  this  research  is  to  make  a  detailed  investigation  of 
finite-amplitude  standing  waves  in  air  at  ambient  tempera- 
tures in  cavity  configurations  producing  degenerate  or 
nearly  degenerate  modes,  and  compare  the  results  to  the 
present  state  of  the  theory  of  Coppens  and  Sanders. 
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II .    BACKGROUND  AND  THEORY 

The  study  of  distortion  of  intense  acoustic  waves  began 
in  1868  when   Kirchoff  [5]  used  the  nonviscous  hydrodynamic 
equations  to  predict  a  change  in  waveform  for  a  plane  travel- 
ing longitudinal  wave.   The  investigations  of  finite  ampli- 
tude effects  in  traveling  waves  continued  [6,7,8],  and  in 
1954  Keller  [9]  was  able  to  extract  solutions  for  standing 
waves  in  a  closed  tube.   His  investigations,  based  on  a  non- 
dissipative  medium,  yielded  useful  results  for  frequencies 
away  from  resonance,  but  predicted  infinite  amplitudes  at 
resonance . 

In  absorptive  media,  only  those  waves  of  relatively  high 
amplitude  will  experience  a  change  in  form.   The  solution  of 
this  nonlinear  wave  equation  requires  second  order  terms  or 
higher,  which  Keck  and  Beyer  [10]  developed  using  perturba- 
tion methods  and  Fourier  series  representation  of  the  wave- 
form.  Coppens  and  Sanders  [1]  extended  the  perturbation 
approach  to  include  wall  losses  predicted  by  Ray leigh-Kirchof f 
[11] ,  and  showed  excellent  agreement  with  experiments  con- 
ducted in  a  rigid-walled  tube  at  low  levels  of  nonlinear 
interactions.   The  experimental  work  by  Beech  [12]  with 
supporting  calculation  by  Ruff  [13]  showed  that  at  high 
excitation  levels,  a  significant  difference  developed  between 
theory  and  experiment.   Winn  [14]  experimentally  demonstrated 
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that  the  Ray leigh-Kir chof f  loss  mechanisms  were  not  com- 
pletely suitable  in  describing  the  phase  relationships  ob- 
served in  the  distorted  waveform.   Coppens  and  Sanders  [3] 
then  revised  their  model  to  utilize  empirically  determined 
losses  and  resonant  frequencies,  thereby  extending  the 
excellent  agreement  between  theory  and  experiment  to  higher 
excitation  levels.   A  further  extension  of  their  model  to 
two-  and  three-dimensional  standing  waves  in  cavities  was 
purposed  and  experimentally  found  by  Lane  [3] ,  to  success- 
fully predict  the  major  features  of  the  harmonic  content  for 
finite-amplitude  standing  waves  in  rectangular  cavities. 
Devall  [4]  ,  however,  found  that  if  degeneracies  existed, 
the  model  failed  to  account  for  the  experimentally  observed 
excitation  of  non-family  modes. 

The  model  of  Coppens  and  Sanders  is  based  on  a  three- 
dimensional,  nonlinear  wave  equation  with  a  dissipative 
term  describing  absorptive  losses  encountered  by  plane  stand- 
ing waves  in  rigid  walled  cavities.   In  one  dimension,  the 
wave  equation  takes  the  form 


00  g 

1  (- 


2       2     2 
n=l  9x     C  Z     9t 
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-  2  9^         _ 

+   D   )  TT-Z     =  b   d 
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where 
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C   =  particle  displacement  due  to  n    harmonic 
n 


b   =  1/2(1  +  C  /C  ) 
p  v' 
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C   =  specific  heat  for  medium  at  constant  pressure 
P 

C   =  specific  heat  for  medium  at  constant  volume, 
v 

The  term   D    depends  on  the  loss  mechanism , and  for  Rayleigh 

n 

-Kirchoff  wall  losses  is 


D   =  6.  ( 

n     1 


.  3/2    2~. 

(to   )  n      dx  dt 


1/2     2 

n      dx 


where 
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a  coefficient  containing  the  shear  viscosity  and 
thermal  conductivity  of  the  medium 
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frequency  (angular)  of  fundamental  resonance 
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and   co    -co,     is  the  difference  of  the  half  power  frequen- 
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cies  above  and  below  the  fundamental. 

Equation  1  can  be  reformulated  [3]  to  introduce  the  em- 
pirical resonance  frequencies  and  absorption  coefficients. 
This  results  in 
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(U)    -  to    is  the  difference  of  the  half  power  frequencies 
n  + 

above  and  below   u)   .) 

nr 


-7T-2-  =  2-—  =  M    T,  R   cos  K  (L-X)  cos  n(u)t  +  6  ) 

ox       _  2  n        n  n 
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p   =  acoustic  pressure 

p   =  density  of  medium 
o 

C   =  thermodynamic  speed  of  sound  in  medium 
o 

2 

M   =  mach  number  =  P,/p  C 

loo 

P   =  peak  pressure  of  fundamental  component  of  wave 

r-     r-     ■  ■  r-  til    , 

R   =  Fourier  coefficient  of  n    harmonic  component, 
n  . 

normalized  such  that  R   =  1 


and 


<f)   =  phase  angle  of  n    harmonic  component,  where  (j>  =  0. 
n  1 

Equation  2  can  be  changed  into  a  time- independent  form,  which 
when  solved  for  boundary  conditions  of  rigid  walls  produces 

,(,  n=l 
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A  quantity  E(n)  is  defined 


E(n)  = 


0)    -  no), 
nr       lr 


nw 


lr 


to  indicate  the  position  of   0)     relative  to  the  ideal  har- 

nr 

monic  frequencies,  no)    .   Values  for  Q(n)  and  E(n)  are 
directly  determined  from  infinitesimal-amplitude  analysis  of 
the  cavity,  and  are  used  as  computer  input  parameters. 
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Ill .    EXPERIMENTAL  CONSIDERATIONS 

A.    APPARATUS 

The  rectangular  cavity  used  in  this  research  (Figure  1) 
was  constructed  from  0.982  in.  milled  aluminum  plates.   The 
interior  of  the  cavity  was  12.002  in»  long  and  2.502  in.  high, 
while  the  width  could  be  varied  from  approximately  5.50  in. 
to  7„00  in.  in  0.25  in.  increments.   All  joints  were  right 
angles,  to  which  a  thin  layer  of  silicon  grease  was  applied 
prior  to  assembly  to  hermetically  seal  the  cavity.   Table  I 
shows  the  theoretically  predicted  normal  mode  frequencies 
for  ideal  rectangular  cavities  corresponding  to  the  seven 
nominal  configurations. 

Acoustic  waves  were  introduced  into  the  cavity  by  means 
of  a  piston  located  in  a  2.25  in.  diameter  port  in  the  floor. 
A  single  lubricated  0-ring  was  used  to  produce  a  seal  between 
the  piston  and  the  porto 

On  a  side  wall  opposite  the  piston,  another  port  was  cut 
for  a  1/4  in.  diameter  Bruel  and  Kjaer  type  4136  condenser 
microphone.   This  port  contained  two  small  O-rings  to  provide 
a  seal  for  the  microphone.   The  axes  of  the  piston  and  micro- 
phone were  mutually  perpendicular  to  minimize  coupling  of 
mechanical  vibration. 

A  block  diagram  of  the  apparatus  is  shown  in  Figure  2. 
The  piston  was  driven  by  an  M-B  Electronics  Model  EA1500 
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Exciter  which,  in  turn  was  driven  by  two  M-B  Electronics 
Model  2120MB  Power  Amplifiers  operating  in  parallel.   The 
driving  signal  was  produced  by  a  General  Radio  1161-A  Coher- 
ent Decade  Frequency  Synthesizer,  configured  to  provide  fre- 
quencies precise  to  0.001  Hz.   A  General  Radio  Model  1160P2 
Sweep  and  Marker  Generator  and  Model  GR1569  Automatic  Level 
Regulator  made  automatic  collection  of  data  possible. 

The  piston  movement  was  continuously  monitored  by  means 
of  an  Endevco  Model  2215  Accelerometer  mounted  within  the 
piston.   The  output  voltage  of  this  accelerometer  was  mea- 
sured on  a  Hewlett  Packard  Vacuum  Tube  Voltmeter  Model  HP400D, 
and  the  output  of  this  meter  displayed  on  a  Fairchild  Model 
776H  Dual-trace  Oscilloscope.    (The  input  impedance  of  the 
oscilloscope  was  too  low  to  allow  the  accelerometer  to  be 
directly  connected  to  the  oscilloscope.) 

The  sound  pressure  in  the  cavity  was  sensed  by  the  Bruel 
and  Kjaer  type  4136  condenser  microphone  with  matching  pre- 
amplifier type  2801.   The  output  of  this  preamplifier  went  to 
five  other  pieces  of  equipment.    (1)  A  Hewlett  Packard  HP400D 
VTVM  to  measure  the  overall  voltage  level.  (2)  Second  input 
of  the  Fairchild  Oscilloscope.   This  allowed  a  continuous 
display  of  the  acoustic  signal  present  in  the  cavity.   The 
oscilloscope  was  externally  triggered  by  the  output  of  the 
synthesizer.   (3-5)  Three  Hewlett  Packard  HP302A  Wave  Analyzers 
with  7Hz  bandwidth  set  to  AFC  (Automatic  Frequency  Control) , 
so  that  they  automatically  followed  the  selected  signal  as 
the  frequency  slowly  varied. 
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Two  of  the  wave  analyzers  were  set  to  the  second  and 
third  harmonic  of  the  sound  signal.   Their  outputs  were 
passed  to  two  Clevite  Brush  Recorders  Model  220,  and  a  time 
record  of  the  second  and  third  harmonic  signal  contents 
were  printed  on  one  channel  of  the  paper  tape  of  their  res- 
pective recorders. 

The  output  of  the  third  wave  analyzer,  set  to  the  funda- 
mental frequency,  was  used  as  a  control  signal  for  the  Auto- 
matic Level  Regulator.   The  Regulator  controlled  the  driving 
level  to  the  2120MB  Power  Amplifiers,  so  that  the  level  of 
the  first  harmonic  remained  constant  while  the  frequency 
was  swept  through  a  cavity  resonance.   The  Sweep  and  Marker 
Generator  was  used  to  automatically  sweep  the  frequency  about 
a  set  center  frequency.   The  Sweep  Generator  has  various 
sweep  speeds  and  produces  a  marker  signal  at  the  center  fre- 
quency,as  well  as  at  set  intervals.  The  marker  signal,  ampli- 
fied by  a  HP400D  VTVM,  was  passed  to  the  second  channels  of 
both  Brush  Recorders.    (It  was  necessary  to  amplify  the 
marker  signal  because  the  Sweep  Generator  was  operated  at 
its  lowest  speed,  where  the  marker  voltage  was  insufficient 
to  produce  a  mark.) 

The  final  element  in  the  experimental  setup  was  a  control 
and  safety  circuit  specifically  designed  and  constructed  for 
this  project.   This  circuit  performed  three  basic  tasks. 
The  first  was  a  switching  function  that  allowed  connection 
of  all  the  appropriate  signals  to  the  equipment  so  the 
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Automatic  Level  Regulator  maintained  control  of  the  drive 
level.   The  opposite  switch  position  bypassed  the  Regulator 
completely  and  allowed  the  drive  signal  to  be  controlled  at 
the  Frequency  Synthesizer. 

The  second  function  of  this  circuit  was  to  provide  a 
safety  relay.   Since  the  output  of  the  Level  Regulator  is 
inversely  proportional  to  the  voltage  of  the  control  signal, 
loss  of  that  signal  causes  the  Regulator  to  output  its  maxi- 
mum voltage  to  the  2120MB  Power  Amplifiers.   To  prevent 
possible  damage,  the  safety  circuit  monitors  the  signal  to 
the  power  amplifiers  and  removes  this  signal  wherever  it 
exceeds  a  value  set  by  the  Associated  Hewlett  Packard  HP467A 
Power  Amplifier  as  shown  in  Figure  2.   Once  the  circuit  is 
tripped,  indicated  by  a  panel  light,  it  remains  open  until 
reset  by  the  safety  circuit/bypass  switch. 

The  last  job  performed  by  the  control  and  safety  circuit 
is  that  of  a  variable  attenuator,  which  allows  the  operating 
voltage  levels  of  the  Regulator  to  be  maintained  near  its 
mid-scale  range.   This  also  provides  an  inherent  safety 
feature  by  limiting  the  maximum  drive  signal.   This  feature 
could  also  be  bypassed. 

B.    STRENGTH  PARAMETER  AND  MICROPHONE  SENSITIVITY 
The  strength  parameter  defined  as 

SP  =  MbQx  , 
where 

M  =  P1/pc2 
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is  the  basic  quantity  which  characterizes  the  strength  of 
the  finite-amplitude  interaction.   To  determine  SP,  it  is 
necessary  to  know  P  ,  which  in  turn  can  be  calculated  from 
the  microphone  output  voltage  if  the  microphone  sensitivity 
is  known.   Microphone  sensitivity,  obtained  with  a  Breul  and 
Kjaer  Model  4220  pistonphone,  was  found  to  be 


Sw  =  (1.56  +  0.03)  x  10~3  VOLT/ (nT/M2) . 
M  — 


Using  the  relationships 


and 


S   =  V  /P 

M     M   rms 


P,  =  /2  P 

1        rms 


one  can  rewrite  the  strength  parameter  as 


-3 
SP  =  (7.62  x  10   )  V  On  , 

M  1 


where   V    is  the  rms  voltage  of  the  fundamental 
M 


C.    FREQUENCY  PARAMETER 

The  frequency  parameter,  defined  by 


FP  =  2(fn  -  nf1)  Q1/nf]_  , 


normalizes  the  driving  frequency  to  the  fundamental  resonance 
of  the  system.   For  example,  FP  equal  to  +  1.00  corresponds 
to  driving  the  cavity  at  the  1/2-power  point  of  the  funda- 
mental resonance. 
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D.  AUTOMATIC  DATA  COLLECTION 

Previous  investigations  [3,4,12,14]  encountered  a  drift 
of  the  resonant  frequencies.   Slight  temperature  changes 
were  postulated  to  be  the  cause.   To  minimize  this  effect 
the  system  was  allowed  to  warm  up  by  driving  the  cavity  for 
four  hours  at  a  strength  parameter  approximately  equal  to 
those  to  be  investigated.   This  significantly  reduced  the 
rate  of  frequency  drift.   This,  coupled  with  the  10-  to  20- 
fold  decrease  in  run  time  that  automatic  data  collection  pro- 
vided, greatly  reduced  frequency  drift  from  0.8Hz  per  run 
to  less  than  0o05Hz  per  run. 

Another  advantage  of  this  automated  system  was  the 
simultaneous  recording  of  harmonics  as  a  continuous  function 
of  frequency.   The  number  of  harmonics  investigated  during 
any  one  run  was  limited  only  to  the  number  of  available 
HP302A  Wave  Analyzers  and  Brush  Recorder  channels. 

E.  SYSTEM  ALIGNMENT 

The  primary  factor  in  setting  the  lower  limit  to  the 
measurement  of  harmonic  content  was  the  amount  of  distortion 
injected  into  the  cavity  by  the  piston,,   Misalignment  of  the 
piston  within  its  port  causes  distortion  of  the  piston  motion 
which  was  detected  by  analyzing  the  harmonic  content  of  the 
accelerometer  output  voltage.    (With  the  piston  removed  from 
its  port,  distortion  levels  were  of  the  same  order  of  magni- 
tude as  the  general  electronic  noise  level  of  the  accelero- 
meter system  when  at  rest.)   Before  each  run,  the  piston  was 
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aligned  within  its  port  so  that  the  acceleromete r  output 
voltage  contained  less  than  0.1%  second  harmonic  and  0.2% 
third  harmonic. 

After  testing  various  configurations,  the  one  selected 
for  ease  of  alignment  consisted  of  the  cavity  mounted  on  a 
tripod  support  directly  above  the  exciter,  which  was 
supported  by  its  own  base  in  an  upright  position.   Both  the 
cavity  tripod  and  the  exciter  rested  on  a  5/16  in.  pad  of 
stiff  acoustic  rubber.   Once  the  system  was  aligned,  only 
an  occasional  minor  adjustment  was  necessary  to  maintain 
low  piston  distortion. 
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IV.    DATA  COLLECTION  PROCEDURES 

A.    PRERUN  PROCEDURES 

Because  the  system  was  designed  to  collect  the  bulk  of 
the  data  automatically,  careful  attention  to  prerun  proce- 
dures was  essential  for  an  accurate  data  run.   The  system 
was  warmed  up  for  at  least  four  hours  prior  to  any  final  run 
preparations,  so  as  to  keep  any  frequency  drift  to  a  mini- 
mum.  After  this  warm  up  period  and  prior  to  anything  else, 
the  piston  motion  was  analyzed  using  the  HP302A  Wave 
Analyzers,  and  the  system  aligned  as  necessary. 

The  HP302A's  were  then  reconnected  to  monitor  the  micro- 
phone output  and  set  to  track  (AFC  position  of  function 
switch)  the  fundamental,  second,  and  third  harmonic  frequen- 
cies.  Next,  the  Frequency  Synthesizer  was  set  to  the 
resonance  frequency  of  the  gravest  member  of  the  family  to 
be  investigated,  and  the  signal  amplitude  was  adjusted 
until  the  value  of  the  fundamental  voltage  reached  the 
desired  level.   The  ratios  of  each  harmonic  to  the  funda- 
mental was  calculated  from  the  wave  analyzer,  and  an  appro- 
priate scale  was  selected  for  each  Brush  Recorder  so  that 
the  stylus  was  approximately  at  its  maximum  position.  Once 
the  Brush  Recorders  were  set  they  did  not  need  to  be  re- 
calibrated as  long  as  the  fundamental  voltage  was  maintained 
at  the  prescribed  level. 
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The  Sweep  and  Marker  Generator  was  set  to  sweep  a  fre- 
quency band  equivalent  to  +  3.0  in  frequency  parameter  cen- 
tered at  the  resonant  frequency,  and  an  appropriate  marker 
spacing  was  set.   The  sweep  was  momentarilly  stopped  while 
the  Automatic  Level  Regulator  was  switched  into  the  circuit 
and  brought  up  to  the  prescribed  voltage  level.    (Maximum 
stable  regulation  rate  was  10  dB/sec.)   The  automatic  sweep 
was  then  allowed  to  continue,  and  the  output  of  the  2120MB 
Power  Amplifiers  was  carefully  monitored  to  prevent  acci- 
dental overdriving  of  the  EA1500  Excitor.   Sweep  speed  of 
sixty  seconds,  and  tape  speed  of  5mm/sec  were  used  and  the 
system  was  allowed  to  make  two  full  sweeps,  at  which  time 
Recorder  calibration  and  stability  of  the  fundamental 
voltage  level  was  checked.   If  all  systems  were  properly 
operating,  the  time  was  recorded  and  the  system  made  its 
final  two  frequency  sweeps  while  the  fundamental  voltage 
level  was  continuously  monitored  for  any  undesired  change. 

B.    RUN  SEQUENCE 

The  run  officially  began  when  the  final  two  frequency 
sweeps  commenced.   (All  previous  recordings  were  discarded 
to  prevent  confusion  due  to  possible  frequency  drift.)   Upon 
completion  of  the  final  two  frequency  sweeps,  the  second 
phase  of  the  measurements  followed  immediately.   The  Auto- 
matic Level  Regulator  and  Sweep  Generator  were  switched  out 
of  the  system,  and  the  Frequency  Synthesizer  was  used  alone 
to  carry  out  the  infinitesimal  amplitude  analysis.   Detailed 
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measurements  of  the  Q  and  E  for  the  fundamental,  each  re- 
corded mode,  and,  where  possible,  modes  of   higher  order 
were  performed.   In  measuring  the  Q  ,  it  was  necessary  to 
maximize  the  microphone  output  voltage  by  varying  the  fre- 
quency.  Then  the  frequency  was  shifted  off  resonance  until 
the  voltage  was  -3dB  below  the  maximum.   The  difference 
between  the  -3dB  points  on  both  sides  of  the  maximum  deter- 
mined the  bandwidth  (Af),  and,  assuming  the  resonance  fre- 
quency was  exactly  between  these  two ,  allowed  the  calculation 
of  Q.   In  those  cases  where  degenerate  modes  appeared  and 
-3dB  down  frequencies  could  not  be  achieved,  the  -ldB  and 
-2dB  points  were  measured  for  both  the  family  and  degenerate 
modes  and  their  center  frequencies  and  Q's  could  be  extra- 
polated using  a  method  described  by  Devall  [5] .   Since  Q 
measurements  were  observed  to  be  very  similar  between  the 
runs  for  both  cavity  configurations,  the  average  values  for 
Q  of  the  non-degenerate  modes  were  used  as  the  computer  in- 
put parameters  for  those  modes  where  degeneracies  impaired 
the  measurements.   An  average  set  of  Q  and  E  measurements 
took,  from  five  to  six  minutes.   With  the  present  automatic 
collection  system,  increasing  the  number  of  harmonics  re- 
corded by  the  addition  of  more  analyzers  and  recording 
channels  will  not  effect  the  actual  run  time. 
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V.    RESULTS 

Tables  II  through  V  display  the  properties  of  the 

cavity  determined  from  inf ini tismal-ampli tude  measurements 

for  each  run.   The  tables  represent  one-  and  two-dimensional 

modes  (nOO  and  nnO)  for  the  non -degenerate  and  degenerate 

configurations  (A  and  B  respectively) .   Each  table  lists 

the  strength  parameter  (SP)  of  the  various  runs,  and  contains 

the  center  frequency  of  the  n    harmonic  (f   )  along  with 

On 

their  measured  bandwidths  (Af),  Q's,  and  E's. 

The  theoretical  values  for  harmonic  content  were  calcu- 
lated using  the  values  found  in  these  tables.   Cavity  con- 
figuration B's  design  specifications  required  the  200  and  010, 
400  and  020,  500  and  320,  and  220  and  410  modes  to  be  de- 
generate pairs.   Data  were  taken,  and  theoretical  curves 
generated  in  both  configurations  for  strength  parameters  of 
about  0.29,  0.38,  0.47  (B  only),  0.55,  and  0.65  driving  the 
100  mode  and  0.24,  0.37,  and  0.49  for  driving  the  110  mode. 
In  all  runs  the  piston  was  positioned  with  its  face  flush 
with  the  interior  floor  of  the  cavity. 

Figures  3  through  17  display  the  experimentally  measured 
harmonic  distortion  (dashed  line)  and  the  theoretically  pre- 
dicted values  (dotted  line).   With  few  exceptions,  theore- 
tical and  experimental  values  lie  well  within  experimental 
error  limits.   For  Figures  7,  8,  9,  10,  11,  15,  16  and  17, 
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TABLE    II 


nOO  Modes 

CONFIGURATION 

A  (12"x5, 

(75"x2.5' 

") 

SP 

RUN# 

n 

<H§y 

Af 

Q(n) 

E(nX 
(xl0J) 

(Hz) 

0.295 

5 

1 

568.01 

2.20 

257 

0.0 

2 

1136.64 

3.08 

388 

0.536 

3 

1705.14 

3.87 

439 

0.648 

4 

2273.89 

4.38 

516 

0.809 

(fig. 

,3) 

0.380 

1 

1 

565.64 

2.26 

249 

0.0 

2 

1131.90 

3.08 

367 

0.551 

3 

1697.77 

3.80 

447 

0.507 

4 

2264.17 

4.53 

511 

0.644 

(fig. 

,4) 

0.551 

3 

1 

568.17 

2.36 

241 

0.0 

2 

1137.12 

3.12 

364 

0.688 

3 

1705.76 

3.91 

436 

0.736 

4 

2274.86 

4.40 

517 

0.966 

(fig. 

.5) 

0.651 

4 

1 

568.14 

2.34 

244 

0.0 

2 

1137.03 

3.09 

367 

0.637 

3 

1705.58 

3.87 

440 

0.362 

4 

2274.62 

4.42 

514 

0.912 

(fig. 

.6) 

30 


TABLE  III 


nOO  MODES         CONFIGURATION  B  (12"x6 . 0"x2 . 5" ) 

SP       RUN#    n        f         Af     Q(n)     E  (nX 

(HzT       (Hz)  (x!0J) 


0. 

,292 

9 

1 

566.81 

2.22 

255 

0.0 

2 

, 1134.84 

3.40 

334 

1.08 

3 

1701.66 

3.77 

451 

0.725 

4 

2270.07 

4.43 

512 

1.25 

(fig, 

.7) 

0. 

,381 

11 

1 

566.82 

2.27 

249 

0.0 

2 

1134.90 

3.48 

326 

1.10 

3 

1701.78 

3.86 

440 

0.771 

4 

2270.26 

4.45 

510 

1.31 

(fig. 

.8) 

0, 

.473 

12 

1 

566.94 

2.28 

248 

0.0 

2 

1135.17 

3.47 

326 

1.14 

3 

1702.19 

3.86 

441 

0.808 

4 

2270.83 

4.46 

509 

1.36 

(fig. 

.9) 

0, 

,559 

13 

1 

564.50 

2.31 

244 

0.0 

2 

1130.34 

3.35 

337 

1.19 

3 

1695.06 

3.80 

446 

0.927 

4 

2261.44 

4.39 

513 

1.53 

(fig. 

.10) 

0, 

,651 

14 

1 

564.93 

2.32 

244 

0.0 

2 

1131.23 

3.45 

328 

1.22 

3 

1696.30 

3.86 

439 

0.887 

4 

2263.05 

4.45 

508 

1.47 

(fig, 

.11) 

—  •- 

— 
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TABLE  -IV 

nnO  MODES         CONFIGURATION  A  (12"x5 . 75"x2 . 5") 

SP      RUN#   n-      fn        A£  Q(n)     E(nX 

(Hz)  (xl(T) 


(Hg? 


0.255      8 

1 

1309.28 

3.91 

335 

0.0 

2 

2629.16 

5.26 

486 

0.608 

3 

3930.70 

7.43 

519 

0.725 

(fig. 12) 

0.382      6 

1 

1310.18 

3.92 

334 

0.0 

2 

2621.67 

5.58 

481 

0.499 

3 

3932.90 

7.63 

516 

0.599 

(fig. 11) 

0.504      7 

1 

1308.12 

3.96 

330 

0.0 

2 

2617.87 

5.29 

494 

0.620 

3 

3927.48 

7.45 

527 

0.793 

(fig. 12) 
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TABLE  V 


nnO  MODES         CONFIGURATION  B  (12"x6 . 0"x2 . 5" ) 

SP       RUN#    n        fn        Af     Q(n)     E  (nX 

(Hz?       (Hz)  (xlOJ) 


0.240     15 
(fig. 15) 

1 
2 
3 

1263.38 
2527.79 
3793.84 

4.01 

5.20* 

7.04 

315 

486* 

538 

0.0 

0.406 

0.975 

0.360     16 
(fig. 16) 

1 
2 
3 

1265.96 
2532.95 
3801.56 

4.01 

5.21* 

7.05 

315 

486* 

539 

0.0 

0.404 
0.967 

0.476     18 
(fig. 17) 

1 
2 
3 

1267.64 
2536.31 
3803.72 

4.11 

5.22* 

6.78 

307 

486* 

552 

0.0 

0.403 

0.207 

*These  values  were  obtained  from  averaging  the  Q's 
for  the  same  mode  from  configuration  A. 
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arrows  indicate  the  locations  of  the  associated  degenera- 
cies.  Tables  VI  and  VII  list  the  resonance  frequencies 
and  Q's  of  the  driving  mode,  the  center  frequencies  and  the 
value  of  FP  at  which  the  mode  would  peak  for  each  degenerate 
pair.    (A  perfectly  tuned  cavity  with  no  wall  losses  would 
exhibit  a  frequency  parameter  of  zero  for  these  frequencies.) 
One  point  of  considerable  interest  is  that  on  all  figures 
for  the  degenerate  configuration  (Figures  7-11,  15-17), 
there  is  no  observable  deviation  of  the  experimental  and 
theoretical  curves  for  the  second  harmonic.   This  is  true 
for  the  one-dimensiona lly  driven  mode  (100)  as  well  as  the 
two-dimens iona lly  driven  one  (110). 

The  only  significant  discrepancies  appear  in  Figures  5, 
6,  9,  and  10,  near  the  peaks  of  the  curves  for  the  third 
harmonics.   The  theory  predicts  a  significantly  flat  res- 
ponse for  these  higher  strength  parameters  that  was  not 
observed  experimentally.   It  is  also  interesting  to  note 
the  peak  position  of  the  fifth  harmonic  and  its  effect  on 
the  fourth. 
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VI.    CONCLUSIONS 

The  experimental  system  used  in  testing  the  theory  of 
Coppens  and  Sanders  has  proven  to  be  accurate,  efficient, 
and  very  reliable.   System  limitations  are  few;  however, 
possible  equipment  improvements  are  outlined  in  Appendix  B. 
The  system  is  easily  adaptable  to  cavities  or  tubes  of 
almost  any  shape. 

It  is  significant  to  note  that  there  were  no  observable 
differences  between  the  theory  of  Coppens  and  Sanders  as 
modified  to  include  empirical  losses,  for  either  degenerate 
or  non-degenerate  configurations.   The  theory  successfully 
predicts  the  basic  harmonic  content  for  f ini te -ampl i tude 
standing  in  a  non-perfect  rectangular  rigid-walled  cavity 
when  the  piston  and  its  port  are  relatively  small  perturba- 
tions in  the  cavity's  geometry,  regardless  of  the  existences 
of  degeneracies. 

The  only  deviation  between  theory  and  experiment  can 
be  seen  best  in  Figures  5,  6  and  10.   The  predicted  harmonic 
distortion  for  the  third  harmonic  in  these  cases  shows  an 
apparent  yielding  of  power  to  the  fourth  harmonic  as  strength 
parameter  is  increased.   The  experimental  data  do  not  con- 
form this  feature.   It  is  unclear  as  to  the  reason  for  this 
discrepancy,  but  appears  to  be  an  artifact  brought  about  by 
the  way  the  computer  program  produces  E(n)  for  n's  greater 
than  those  actually  measured.   The  program  presently 
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calculates  values  for  E(n)  not  empirically  determined,  and 
assigns  them  alternating  positive  and  negative  values  of 
magnitudes  consonant  to  those  actually  measured.   In  all 

cases  investigated,  the  experimental  values  for  E(n)  were 

-4 
positive  and  of  order  10    for  the  one  dimensional  cases  in 

question.   It  is  probable  that  a  slight  change  in  the  method 

of  calculating  higher  order  values  of  E(n)  could  alleviate 

this  lack  of  agreement.   A  quick  computer  run  was  made 

using  the  values  obtained  for  run  number  4,  and  a  fabricated 

value  of  E(5)  which  was  positive  and  had  a  magnitude  equal 

to  the  average  of  E(2),  E(3),  and  E(4),  and  it  was  found 

the  flattening  in  the  predicted  values  vanished  and  more 

closely  resembled  the  experimental  data. 
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APPENDIX  A 
INVESTIGATION  OF  TRAPEZOIDAL  CAVITY 

One  feature  incorporated  into  the  design  of  the  rec- 
tangular cavity  was  the  ability  to  turn  the  adjustable  wall 
around.   The  side  of  the  wall  that  was  out  when  in  the  rec- 
tangular configuration,  was  milled  on  an  angle  that  had 
0.122  in.  rise  over  the  12.002  in  length  (Figure  Al).   The 
mean  interior  width  was  adjustable  from  approximately 
5.50  in.  to  7.00  in.  in  0.25  in.  increments.   One  configura- 
tion was  investigated  in  the  same  manner  as  were  the 
rectangular  ones;  a  mean  width  of  5.740  in.  was  utilizedc 

The  cavity  resonances  were  found  to  correspond  closely 
to  the  theoretically  predicted  frequencies  for  an  ideal 
rectangular  cavity  with  an  interior  width  dimension  of 
5.75  in,  but  always  were  displaced  toward  higher  frequencies 
The  cavity  was  excited  at  its  lowest  resonance  frequency 
and  the  harmonic  content  recorded.   Table  A-l  lists  the 
strength  parameters,  center  frequencies,  Q (n ) 's  and  E (n) 's 
measured  in  this  configuration.   It  is  interesting  to  note 
that  these  Q  (n)   *s  were  generally  lower  than  those  found  for 
configuration  A  (Table  II). 

Theoretically  predicted  curves  were  generated  using 
the  values  in  Table  A-l,  and  compared  to  the  observed  har- 
monic content  (Figures  A2-A4).   It  is  clear  that  there  is 
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TABLE    A-l 


nOO  MODES 

TRAPEZOIDAL  CONFIGURATION 

SP      RUN# 

n 

f0n 

Af 

Q(n) 

E(n) 

(xlO3) 

(Hz) 

(Hz) 

0.269     21 

1 

565.67 

2.40 

235 

0.0 

2 

1132.01 

3.18 

355 

0.684 

3 

1698.09 

3.84 

422 

0.636 

4 

2264.69 

4.64 

488 

0.886 

(fig.A2) 

0.456     20 

1  ' 

565.51 

2.36 

239 

0.0 

2 

1131.75 

3.14 

360 

0.642 

3 

1697.68 

4.01 

423 

0.677 

4 

2264.24 

4.61 

491 

0.969 

(fig. A3) 

0.549     19 

1 

565.35 

2.36 

240 

0.0 

2 

1131.47 

3.15 

359 

0.681 

3 

1697.27 

4.05 

419 

0.716 

4 

2263.72 

4.64 

488 

1.02 

(fig.A4) 

55 


EXPERIMENTAL 
THEORETICAL 


nOO  MODE 
TRAPEZOID 
SP=  0.269 
RUN#  21 


0.1 


/ 


X 


\ 


\ 


N 


/: 


\ 


0.05 


/■ 


X. 


/< 


/• 


x. 


S 


/. 


X. 


.s- 


n=2       ,<v 


/.' 


s.- 


Vv 


/ 


o.oi. 


/ 


/ 


V. 
NT- 


S' 


V- 


.•/ 


0.005 


1 


•I 


v. 

V. 

\\ 

.    v. 


:/ 


\- 


.7 


\ 


n=3 


./ 


\ 


y 


y. 


* 


/ 


/ 


y 


w 


n=4 


0.001 


i t^i 


n=5 


+  1 


+  2 


FREQUENCY    PARAMETER    (FP) 
FIGURE   A2. 


+  3 


56 


0.1 


0.05 




EXPERIMENTAL 

THEORETICAL   , 

y. 

-  " 

/  .' 

y  y 

y .' 

/  y 

yy 

y  y 

/  .• 

n=2 

y  y 

S 

#  • 

o.oi. 


0.005 


\ 


nOO  MODE 
TRAPEZOID 
SP=  0.456 
RUN#  20 


\ 


v. 


\ 


\ 


V 


X 


\ 


X 


-      / 


/"~*s 


S 


N 


/ 


/.' 


/: 


i: 


/: 


/• 


/• 


n=3 


/• 


/. 


/• 


«    • 


V 


V 


V 


V. 


\. 

v 


V 


V. 


K. 


\. 


\. 


\ 


\ 


\ 


*.\ 


n=4 


0.001 


n=5 


i i 


0  +1 

FREQUENCY   PARAMETER    (FP) 
FIGURE   A3. 


+  2 


+  3 


57 


0.1 


0.05 


n 


EXPERIMENTAL 

THEORETICAL  /.' 

/.' 
/.' 

/.* 

A' 

sy 

/.- 

n=2  y> 

y 


y. 


N 


V 


H. 


nOO  MODE 
TRAPEZOID 
SP=  0.549 
RUN#  19 


\. 


V 


S 


\ 


.'/ 


•y 


.\ 


^ 


y 


/ 


7 


/ 


/ 


/ 


/ 


/ 


/.' 


/: 


/. 


0  .  0  1-  -    n=3 


A" 


0.005 


0.001 


/ 


n=4 


;  n=5 


i i i 


'\ 


V. 


V 


V. 
V 


V. 
V. 


\. 


\ 


* 


\ 


\ 


<•! 


FREQUENCY    PARAMETER    (FP) 
FIGURE   A4. 


+  2 


+  3 


58 


excellent  agreement  in  all  three  runs.   The  same  phenomena 
exhibited  in  Figures  5,  6,  10,  and  11  is  seen  in  Figure  A4, 
and  is  believed  to  be  an  artifact  of  the  computer's  method 
of  calculating  values  of  E(n)  of  higher  order  than  those 
empirically  determined. 

It  appears  that  the  theory  of  Coppens  and  Sanders  can 
successfully  predict  harmonic  content  for  cavities  slightly 
perturbed  from  rectangular  configurations. 
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APPENDIX  B 


SUGGESTIONS 


The  major  limitation  encountered  in  this  research  was 
the  inability  to  record  more  than  2  harmonics  at  one  time. 
What  is  needed  is  a  six- to  eight- track  recorder  and  at 
least  two  additional  Hewlett  Packard  HP302A  Wave  Analyzers. 
All  harmonics  to  be  investigated  must  be  recorded  simul- 
taneously to  keep  any  frequency  drift  from  becoming  a  factor 
in  these  measurements. 

The  General  Radio  Sweep  and  Marker  Generator  model 
1160P2  was  borrowed  from  a  separate  activity.   Use  of  this 
piece  of  equipment  is  essential  to  the  operation  of  this 
system . 

The  final  suggestion,  but  of  no  less  importance,  is  the 
necessity  of  the  0.00  through  0.09  Hz  digital  component  for  the 
existing  General  Radio  model  1161-A  Coherent  Decade  Frequency 
Synthesizer.   The  accuracy  to  which  the  Sweep  Generator's 
center  marker  is  located  is  dependent  on  the  Synthesizer's 
digital  components.   This  final  component  would  eliminate 
any  marker  ambiguity  that  is  introduced  in  the  system's 
present  configuration. 
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